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Abstract 

We consider a class of higher order corrections with arbitrary power n of the 
Q I curvature tensor to the standard gravity action in arbitrary space-time dimension 



D. The corrections are in the form of Euler densities and are unique at each n and 
D. We present a generating functional and an explicit form of the corresponding 
conserved energy-momentum tensors. The case of conformally flat metrics is 
discussed in detail. We show that this class of corrections allows for domain wall 
solutions since, despite the presence of higher powers of the curvature tensor, the 
singularity structure at the wall is of the same type as in the standard gravity. 
However, models with higher order corrections have larger set of domain wall 
solutions and the existence of these solutions no longer depends on the presence of 
cosmological constants. We find for example that the Randall-Sundrum scenario 
can be realized without any need for bulk and/or brane cosmological constant. 
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The idea of calculation of the effective quantum corrections to the energy-mo- 
mentum tensors in the nontrivial gravitational background is relatively old |ll|, 0]. If 
gravity is treated classically such a calculation corresponds to a one-loop result in 
quantum theory. If the background metric is conformally flat then there is no massless 
particle creation and the result for the vacuum polarization is local, therefore the 
use of the effective lagrangian is justified. These corrections are in the form of an 
expansion in powers of the curvature tensor. They belong to two classes. One of them 
("topological") consists of terms that can be obtained from an effective action in the 
first order formalism without need for any metric. The second class ("non-topological" ) 
consists of terms coming from an action that necessarily involves a metric. 

In this paper we are interested in the first class corresponding to the set of Euler 
densities of arbitrary order n {n being a power of the curvature tensor) 0. When 
the dimension of the space-time D is equal to 2n the Euler density gives part of the 
conformal anomaly. This type of anomaly is unique and preserves scale invariance |^, 
(the other type requires introduction of a scale through regularization and vanishes for 
vanishing Weyl tensor). It is important to stress that these higher order gravitational 
terms should always be included in the presence of the quantum matter fields - for 
example in the case of the Standard Model fields[^ 0. 

There is another, independent motivation to consider terms of higher order in the 
curvature tensor. Such terms appear in the a' expansion of the string theory effective 
action and it was shown that the quadratic terms can be put in the form of the Gauss- 
Bonnet combination {n = 2 Euler density) . 

It is interesting to note that Euler densities are the only higher order terms in the 
action that do not introduce ghosts and assure "Cauchy-like" evolution of the initial 
conditions. 

From the form of the Euler densities it follows that for arbitrary D and n all the 
energy-momentum tensors belonging to the first class do not have derivatives acting on 
the curvature tensor. For the case of the conformally fiat metrics we find a generating 
functional for these energy-momentum tensors. We analyse the domain wall solutions 
in this class of models. We show that such solutions can be constructed since, even in 
the presence of higher powers of the curvature tensor, the type of singularity at the wall 
is the same as in the standard gravity. In the usual domain wall metric the curvature 



has a delta-like singularity on the wall. Thus, one would naively expect that the higher 
powers of the curvature tensor introduce meaningless expressions (higher powers of the 
Dirac delta function) which can not correspond to any sources. The Euler densities 
seem to be the only combinations of the higher order terms that do not lead to such 
meaningless expressions and which allow for the domain wall metric with the usual 
sources. 

The number of the domain wall solutions with higher order terms in the action 
is larger than in the standard case and the cosmological constant is in general no 
longer crucial for their existence. It may even happen that the Randall-Sundrum type 
scenario p can be realized without bulk and/or brane cosmological constant. The fact 
that there exist brane-like solutions without any cosmological constant does not solve 
the cosmological constant problem - it is, as usual, the fine tuning problem since the 
solution requires then some relations among the couplings. 

In this paper we assume that all the higher order corrections to the lagrangian are 
of the form of Euler densities which in D dimensions are defined (in the form notation) 
as 

/^"^ = m_^2n)! '°^"^-"^^"'"' ^"' ^"'"""'" ^ e"'"^' A ■ ■ ■ e"^ . (1) 

For D = 2n they are invariants and formally total derivatives. However, careful regu- 
larization does not allow us to discard them neither in the action nor in the equations 
of motion since they correspond to the conformal anomaly. 
If we have a lagrangian of the form 

C = -Y: '^n/^") (2) 

n=0 

then the equations of motion from the variation of the vielbein read (for n < {D — 1)/2) 

y ^4-P~2n) Jia,a2 ^ . ..Jia2„^,a2n ^ a2„+i /\ . . . /\ ^an-i _ Q . (3) 

We can write the curvature two-form as 

R"'' = C''^ + ^^{e''AK^-e^ AK") (4) 

where C"^'' is a two-form composed of the Weyl tensor C^ypfj while K"^ is a one-form 
defined (for invertible vielbeins e^) as K" = Kf^^e^"'dx'' = K"-be^ with 
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In this paper we consider the background metrics for which the Weyl tensor van- 
ishes. In this case the curvature two-form (^) can be expressed in terms of the Ricci 
tensor and the curvature scalar only. 

Multiplying (^ by e^ we get 

= ^ E ^n p r2HD -~2n - 1)! (^"-^"^ • • • ^"-"^° ^ P"^"^-) ^-"'^" • • • ^-" " 

where H^'^^ are defined by this equality and will be described in detail below. 

It is not too difficult to prove that the generating functional for the tensors HJ^^ is 
given by: 

oo 

M^, (t, K) = det( l-tK)(il- tK)-') =Y.e (h^^A (7) 

where from now on we use the matrix notation (1 denotes 5!^ while K denotes K^^). 
Differentiating equality C^„p^ = over x^ and using Bianchi identity we get 

By differentiating M^i, (t, K) over x^ and using (P) it is straightforward to prove that 
tensors H^Jl) are covariantly conserved. A more tedious calculation shows a nontrivial 
fact that (when the Weyl tensor vanishes) H^Jl) are the only covariantly conserved 
tensors which are algebraically composed of the Ricci tensors. 

From the generating functional (0) one gets the recurrence relations for H^J!^: 



HtD = Hj:^-'^K^ - -g,.H^;-'^K^^ (9) 

starting with H^^J = g^^] the exphcit formula (for n > 0) reads 
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where Q = Ylk=i ^1k and the sum is over all sets of n non-negative integers Qk for 
which Q < n. 



As an example we present the first four tensors HJ^^: 



9i^u 



i/W = R --Ra 

-"/ii/ ~ RfipR V — oTjTZri) '^'^ ^>^'^ 1 ~2'^p ~^ Q(r) — I) j ' 

Hj^ = R^^R-pR^^-J^±1-RR^^R\ + R^(-^R^^RP- + £1±^R^rA 

+9t^u [--RraR'^'^Rpr + 4m _ ^N^^P^i?'"" " 48m _ 1)2 ^J " ^^^) 

HJ^J was used (for D = 4) in [0] to prove tliat the Einstein-Hilbert gravity with such 
a correction does not have a cosmological singularity and can describe inflation. 

One can see from the definition (P that in D space-time dimensions forms /'•"^ are 
nonzero only for < n < [D/2] (therefore also Umax < [-D/2] in (|^)). The situation 
is different for the tensors HJ^\ Using the generating functional (0) or the explicit 
formula (jT^) it is straightforward to show that HJ^^ are nonzero in the larger range 



< n < D. The difference arises because the energy-momentum tensors T^"-* derived 
from I^^' are proportional to H^^ with coefficients vanishing for n > D/2: 



;-2)"-in!(D-n-l) 
' (D-2)"(D-2n-l)!' 






The case with D = 2n should be treated with care because it is related to the conformal 
anomaly. 

Let us illustrate this point by considering n = 2. The Euler density in this case 
is the famous Gauss-Bonnet combination: /(^^ = R^/BpaR'^^'"' - ^.R^pR'^'^ + R^. The 
energy-momentum tensor obtained from this density is given by (^): 

T^S = ^^2 (-4i?"^i?„;./3. + 2R^^p^R,^f^ - AR^^R", + 2RR^, - \9^.J^^^) ■ (13) 

Let us note that this energy-momentum tensor is covariantly conserved even with 
nonvanishing Weyl tensor. Relation ([T3|) for the conformally fiat metric gives 

_ 4(C-4)(J-3) 

-^p.u "-2 /jj _ 2)2 V^^ ^ ^ 



The coefficient in the above equahty contains the factor {D — 4) but it is compensated 
by a factor T{D — A) present in K2 due to regularization of the effective action so T^^/^^ 
should be taken into account. 

Let us now use the above formahsm to investigate the possibihty of warped metric 
solutions in theories with higher order gravity terms in the lagrangian (|^). 

In this paper we will consider a model in D-dimensional space-time with a (D — 2)- 
brane (a domain wall). Its action is the sum of the bulk and brane contributions: 



'bulk ~r '-'brane ; 



•S* — 5'buik + S\j 
.Sbuik = -/rf^Xv/=^ 5] «:„/("), (15) 

Sbrane = " / d^'^xJ-g (A + . . .) 



The metric on the brane is given by g^y{x'') = g^i^{x^, y = 0) where y = x^ and from 
now on /i, z/, . . . = 1, . . . ,D — 1 while M, N, . . . = 1, . . . ,D. In the brane action we 
write explicitly only the most important term - the brane cosmological constant. The 
bulk gravitational interactions are described by the sum of terms J^"-* defined in (|ip. 
The first two terms are known from conventional gravity. The one with n = 1 is the 
usual Hilbert-Einstein term, I^^^ = R, and its coefficient, ki, depends on the fields 
normalization. We will use normalization for which ki = 1 (the other frequent choice 
is (2/t;^)^^). The term with n = corresponds to the cosmological constant: I^^^ = 1, 
kq = A. The maximal number of the higher order terms is Umax ^ [-D/2] as previously 
discussed. 



We obtain the equations of motion by differentiating the action (|1^) with respect 
to ^*^^: 

Rmn - -RgMN = -AguN + -Xgif^JM^N^iy) - H ^MTv • (16) 

^ ^ -^ n=2 

We see that TJ,"^ given by (0), for n > 2 can be treated as contributions to the 
energy-momentum tensor coming from higher order interactions present in the action 

Let us now look for the domain wall solutions in D-dimensional space-time which 
are flat from the {D — l)-dimensional point of view. The metric in this case can be 
written in the form 

ds^ = e-^^'^y^'q^^dx^'dx" + dy'^ . (17) 



Using the generating functional (|^ we can find -f^l/jv ^^ ^^is metric: 

2n-2 



J^MN — 9MN 




dy y ^""JD-Ady^ 



Tin) 



Tliis explicit formula shows that all Hmj^ vanish for n > D and that the second 
derivative of f{y) is absent in Hj^jj and appears only linearly in HJ^. This is quite 
surprising for a theory with higher orders of the curvature tensor in the action and, 
as we will see shortly, is crucial for the existence of the domain wall solutions. This 
property follows from the unique structure of the Euler densities (|1|) and the fact that 
the second derivative of f{y) appears only in -R'dd and is absent in K^^,. Although 
higher powers of d"^ f{y)/dy'^ are present in separate terms on the r.h.s. of eq. (0) they 
exactly cancel in the full sum. 

Using tensors (|T8D we find that the equations of motion (|T6|) for the domain wall 
metric (|T^ reduce to just two conditions: 

where the coefficients Pn are equal 

The absence of the third or higher derivatives of f{y) and/or higher powers of the 
second derivative d"^ f {y) / dy"^ in ([T8|) allows for the solution of (|T6D of the same type 
as in the standard gravity: 

f{y) = (y\y\ (22) 

where the parameter a in this case must satisfy two algebraic conditions coming from 

mm-- 

EPnCr^" = -A, (23) 

n=l ^ 

T^^Pncr'"-' = -^^— ^A. (24) 

^^ 2n-r 8 ^ ^ 



One of these equations can be used to determine the warp factor coefficient a in terms 
of the parameters pn and one of the cosmological constants (bulk, A, or brane, A). The 
second equation is then the (fine tuning) condition for the other cosmological constant. 

Recently there has been a lot of interest in domain wall solutions in theories with 
the space-time dimension equal to 5. The Randall-Sundrum scenario corresponds 
to D = 5 and only two terms with ra = 0, 1 {umax = 1) in the lagrangian. In such 
a case the square of the warp factor exponent, a^, is just proportional to the bulk 
cosmological constant A and the brane cosmological constant A must be appropriately 
fine tuned: aX = —A. 

However, the equations (^^4|) have other very interesting solutions when the higher 



order terms are present (i.e. Umax > 2). The square of the warp factor exponent a"^ 
is no longer just linearly proportional to A. It depends also on the other lagrangian 
parameters /€„. In general there are up to rimax possible values of a^ for a given value 
of A. This is true also for A = 0. Thus, it is possible to have nontrivial warped metric 
solutions with vanishing bulk cosmological constant for any rimax > 2. There is one 
such solution for each positive zero of the polynomial P{x) = J2n=T PnX^- It is also 
possible to solve eqs. (p3|j2^) with A = 0. The domain wall is then supported only 
by higher order corrections to the action and (in distinction to the Einstein-Hilbert 
gravity) there is no need for conventional sources on the brane. 

Even more interesting situation emerges when the zero of P{x) satisfies simultane- 



ously eq. (|2^) . Then, there exists such a warp factor a that eq. (|23|) is satisfied for A = 
and eq. ( p^ is satisfied for A = 0. In such a case, there is a domain wall solution which 
is fiat from the {D — l)-dimensional point of view with both cosmological constants 
vanishing. This, however, can happen only for rimax > 2 (i.e. D > 6) and for some 
very specific values of the lagrangian parameters k„. The fine tuning of parameters Kn 
replaces the fine tuning of the brane tension A used in the original Randall-Sundrum 
scenario. Therefore, it is a different formulation of the cosmological constant problem. 
The domain wall solutions with the first correction (Gauss-Bonnet term) were dis- 
cussed in different contexts in the literature. The 5-dimensional Gauss-Bonnet gravity 
with dilaton was discussed in detail in Pj. The brane solutions with vanishing bulk 
cosmological constant (but with non-vanishing brane cosmological constant) were dis- 



cussed in gravity coupled to the dilaton (see for example P, ITU| where the first reference 



does not discuss the Gauss-Bonnet term). The mechanism of "self-tuning" presented 
in these papers was shown later [|ll[] to be another form of fine tuning. The possibility 
of brane type solutions with vanishing cosmological constants was also disscussed in 



T2|. The idea of employing quadratic corrections as boundary counterterms to the 



AdS action was used fr example in [^, T^\. The brane configurations in the context 



of AdS/CFT correspondence were considered for example in refs. [|1J, |T5[ and in the 
cosmological context in |[T6[| . 

In conclusion, we considered a class of models with higher order gravity corrections 
in the form of the Euler densities with arbitrary power n of the curvature tensor in 
arbitrary space-time dimension D. We found a generating functional and an explicit 
form of the corresponding energy-momentum tensors. The domain wall type solu- 
tions were constructed for the considered class of models. Such solutions are possible 
despite the presence of higher powers of the curvature tensor because the singularity 
structure at the wall is of the same type as in the standard gravity. The number of 
possible domain wall solutions increases with the order of included corrections and the 
cosmological constant is no longer crucial for their existence (as it was the case for 
models with n < 2 previously analysed in the literature). We showed for example that 
the Randall-Sundrum type scenario can be realized without the bulk and/or brane 
cosmological constant. 

This work was partially supported by the Polish KBN grants 5 P03B 150 20 and 2 
P03B 052 16. 
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